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1 Introduction 



The production cross sections for all the processes at hadron-collider experiments are con- 
trolled by strong interaction physics and, hence, by its underlying field theory, QCD (see 
recent overviews in Refs. |T|, 0, 0, ffl). Studies of QCD at the Tevatron and the LHC have 
two main purposes ||, ^ |7]]. First, they are important to test the predictions of QCD, to 
measure its fundamental parameters (e.g. the strong coupling as) and to extract quan- 
titative information on its non-perturbative dynamics (e.g. the distribution of partons in 
the proton). Second, they are relevant to a precise estimate of the background to other 
Standard Model processes and to signals of new physics. 

This contribution is not a comprehensive review of QCD at high-energy hadron collid- 
ers. It is based on a selection of the topics presented in my introductory lectures at this 
Workshop. The selection highlights the QCD subjects that were most discussed during the 
Workshop and includes a pedagogical overview of some of the corresponding theoretical 
tools. 

After the introduction of the general theoretical framework, I summarize in Sect. ^| the 
present knowledge on the parton densities and its impact on QCD predictions for hard- 
scattering processes at the Tevatron and the LHC. In Sect. ||, I then discuss some issues 
related to processes that are sensitive to the gluon density and, hence, to its determination. 
Section ^ presents a dictionary of different approaches (fixed-order expansions, resummed 
calculations, parton showers) to perturbative QCD calculations. The dictionary continues 
in Sect. [|, where I review soft-gluon resummation and discuss some recent phenomenological 
applications of threshold resummation to hadron collisions. 

The QCD framework to describe any inclusive hard-scattering process, 



in hadron-hadron collisions is based on perturbation theory and on the factorization the- 
orem of mass singularities. The corresponding cross section is computed by using the 
factorization formula || 



<?{Pl,P2;Q, {•••}) = X) / dXl dX2 f a / h l ( Xl > A*f) fb/hA X 2, /40 0-ab(Xip 1 ,X 2 P2]Q, {■■■}', ^40 



The colliding hadrons hi and h 2 have momenta pi and p 2 , H denotes the triggered 
hard probe (vector bosons, jets, heavy quarks, Higgs bosons, SUSY particles and so on) 
and X stands for any unobserved particle produced by the collision. The typical scale 
Q of the scattering process is set by the invariant mass or the transverse momentum of 
the hard probe, and the notation {. . .} stands for any other relevant scale and kinematic 
variable of the process. For instance, in the case of W production we have Q = M w and 
{...} = {Q±, y, . . .}, where M w , Q ± and y are the mass of the vector boson, its transverse 
momentum and its rapidity, respectively. 

The factorization formula (^) involves the convolution of the partonic cross sections a a t 
(where a,b = q,q,g) and the parton distributions f a /h(x,fi F ) of the colliding hadrons. If 



hip^ + h 2 (p 2 ) ^ H(Q, {...})+ X , 



(1) 




+ O ((A Q cd/Q) p ) 



(2) 
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the hard probe if is a hadron or a photon, the factorization formula has to include an 
additional convolution with the corresponding parton fragmentation function d a m{z, p, F ). 

The term O ((Aq C d/Q) p ) on the right-hand side of Eq. @ generically denotes non- 
perturbative contributions (hadronization effects, multiparton interactions, contributions of 
the soft underlying event, and so on). Provided the hard-scattering process (|l|) is sufficiently 
inclusive^], a a b is computable as a power series expansion in as(Q 2 ) and the non-perturbative 
contributions are (small) power-suppressed corrections (i.e. the power p is positive) as long 
as the hard-scattering scale Q is larger than few hundred MeV, the typical size of the QCD 
scale A QCD . 

The parton densities f a /h(x, fi 2 F ) are phenomenological distributions that describe how 
partons are bounded in the colliding hadrons. Although they are not calculable in QCD 
perturbation theory, the parton densities are universal (process-independent) quantities. 
The scale \i F is a factorization scale introduced in Eq. @ to separate the bound-state effects 
from the perturbative interactions of the partons. The physical cross section a(pi,p 2 ] Q, {. ■ .}) 
does not depend on this arbitrary scale, but parton densities and partonic cross sections sep- 
arately depend on fi F . In particular, higher-order contributions to &ab(xiPi, X2P2', Q, {■ ■ ■}', 
contain corrections of relative order (as{Q 2 ) In Q 2 / fJL F ) n . If [Lf is very different from Q, these 
corrections become large and spoil the reliability of the perturbative expansion. Thus, in 
practical applications of the factorization formula (0), the scale fip is set approximately 
equal to the hard scale Q and variations of /if around this central value are used to estimate 
the uncertainty of the perturbative expansion. 

The lower limit x m \ n of the integrations over the parton momentum fractions X\ and x 2 , 
as well as the values of X\ and x 2 that dominate the convolution integral in Eq. (0), are 
controlled by the kinematics of the hard-scattering process. Typically we have x min ^ Q 2 /S, 
where S = (P1+P2) 2 is the square of the centre-of-mass energy of the collision. If the hard 
probe is a state of invariant mass M and rapidity y, the dominant values of the momentum 
fractions are xi j2 ~ {Me ±v ) / \fS (see Fig. [1]). Thus varying M and y at fixed y/S, we are 
sensitive to partons with different momentum fractions. Increasing the parton densities 
are probed in a kinematic range that extends towards larger values of Q and smaller values 
of x 1>2 . 



2 Parton densities 

The parton densities are an essential ingredient to study hard-scattering collisions. Once 
the partonic cross sections have been perturbatively computed, cross section measurements 
can be used to determine the parton densities. Then, they can in turn be used to predict 
cross sections for other hard-scattering processes. 

The dependence of the parton densities[]/ a (x, fi 2 ) on the momentum fraction x and their 
absolute value at any fixed scale \i are not computable in perturbation theory. However, 

§ More precisely, it has to be defined in an infrared- and collinear-safe manner. 

"In the following the parton densities of the proton f a / p are simply denoted by f a and those of the 
antiproton are obtained by using charge-conjugation invariance, i.e. f a /p = f a / p = fa- 
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LHC parton kinematics 




Figure 1: The (x, Q 2 ) plane of the parton kinematics for the production of a heavy system 
of invariant mass M and rapidity y at LHC, HERA and fixed-target experiments. 



the scale dependence is perturbatively controlled by the DGLAP evolution equation 

dJJ^p = *dz 2 } )/a(x/ 2) (3) 

a In ^ i Jx z 

The kernels P a b{cxs, z ) are the Altarelli-Parisi (AP) splitting functions. As the partonic 
cross sections in Eq. (|2|), the AP splitting functions can be computed as a power series 
expansion in as'- 

P ab (a S} z) = asP%°\z) + alP^ LO \z) + <4p£ nu >\z) + 0(ag) . (4) 

The leading order (LO) and next-to-leading order (NLO) terms P^ (z) and P^ LO \z) in 
the expansion are known ||3_ 0f| . These first two terms are used in most of the QCD studies. 



Having determined f a (x, Qq) at a given input scale fi = Qo, the evolution equation (Q) can 
be used to compute the parton densities at different perturbative scales \x and larger values 
of x. 

The parton densities are determined by performing global fits [11] to data from deep- 



inelastic scattering (DIS), Drell-Yan (DY), prompt-photon and jet production. The method 
consists in parametrizing the parton densities at some input scale Qo and then adjust- 
ing the parameters to fit the data. The parameters are usually constrained by impos- 
ing the positivity of the parton densities (/ a (x,/i 2 ) > 0) and the momentum sum rule 

(Ealo f a (x,/J, 2 ) = 1). 

The present knowledge on the parton densities of the proton is reviewed in Refs. [|7|, |12 



Their typical behaviour is shown in Fig. 0. All densities decrease at large x. At small x the 
valence quark densities vanish and the gluon density dominates. The sea-quark densities 
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Overview of Parton Distribution Functions of the Proton 




Figure 2: Typical x-shape of the parton densities (set CTEQ5M at Q = 5 GeV). 



also increase at small x because they are driven by the strong rise of the gluon density and 
the splitting of gluons in qq pairs. Note that the quark densities are not flavour-symmetric 
either in the valence sector (u v ^ d v ) or in the sea sector (u ^ d). 

In addition to having the best estimate of the parton densities, it is important to 
quantify the corresponding uncertainty. This is a difficult issue. The uncertainty depends 
on the kinematic range in x and Q 2 . Moreover, it cannot be reliably estimated by simply 
comparing the parton densities obtained by different global fits. In fact, a lot of common 
systematic and common assumptions affect the global-fit procedures. Recent attempts 
to obtain parton densities with error bands that take into account correlations between 
experimental errors are described in Refs. |7|, Some important theoretical uncertainties 



that are still to be understood are also discussed in Ref. 0. 

The overall conclusion is that the quark densities^ are reasonably well constrained and 
determined by DIS and DY processes, while the gluon density is certainly more uncertain 
IT], [H]]. At small x (x^lCT 3 ), the gluon density f g is at present constrained by a single 



process, namely DIS at HERA. Thus, large higher-order corrections of the type (as In l/x) n 
could possibly affect the extraction of f g . Assuming that f g is well determined at small x, 
the momentum sum rule reasonably constrains f g at intermediate values of x (x ~ 10 -2 ). 
Jet production at the Tevatron at low to moderate values of the jet transverse energy Ex 
can also be useful in constraining the gluon distribution in the range 0.05^ x ^0.2. At 
large x (x ~ 10 _1 ), the most sensitive process to f g is prompt-photon production. Since, 
at present, prompt-photon data are not well described/predicted by perturbative QCD 
calculations, they cannot be used for a precise determination of f g . Further discussion on 
these points is given in Sect. [| 

The conclusion that the gluon density is not well known can also be drawn by inspection 
(see Fig. 0) of the differences between the most updated analyses performed by the CTEQ 



"Uncertainties on the determination of the quark densities at very high x are discussed in Refs. |l4|, [l5|, 

ra. 
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Comparison of CTEQ5 and MRS98 gluon distributions 
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Figure 3: Comparison between the gluon densities of the CTEQ and MRST groups. 



Collaboration and the MRST group. 

The differences between the MRST gluons and the CTEQ ones are due to the fact 
that the two groups used different data sets. The various gluon densities are very similar 
at small x, because in this region both groups used the HERA data. The MRST group 
includes prompt-photon data in the global fit: these data constrain the gluon directly at 
x ^ f CT 1 and indirectly (by the momentum sum rule) at x ~ f CT 2 . The CTEQ group does 
not use prompt-photon data, but it includes Tevatron data on the one-jet inclusive cross 
section. These data give a good constraint on f g in the region 10~ 2 < x < f0 _1 . 

There are also differences within the MRST and CTEQ sets. The various gluon den- 
sities of the MRST set correspond to different values of the non-perturbative transverse- 
momentum smearing that can be introduced to describe the differences among the prompt- 
photon data that are available at several centre-of-mass energies. The CTEQ5M and 
CTEQ5HJ gluons correspond to different assumptions on the parametrization of the func- 
tional form of f g (x, Qq) at large x; the CTEQ5M set corresponds to the minimum-x 2 
solution of the fit while the CTEQ5HJ set (with a slightly higher x 2 ) provides the best fit 
to the high-Er tail of the CDF and DO jet cross sections. 

This brief illustration shows that the differences in the most recent parton densities are 
mainly due to either inconsistencies between data sets and/or poor theoretical understand- 
ing of them. A more quantitative picture of the dependence on x and Q 2 of the gluon 
density uncertainty is presented in Fig. |j. 

We can see that the DIS and DY data sets weakly constrain f g for f0 _1 . Since the 
AP splitting functions lead to negative scaling violation at large x, when f g (x, Q 2 ) is evolved 
at larger scales Q according to Eq. (Q) the gluon uncertainty is diluted: it propagates at 
smaller values of x and its size is reduced at fixed x. 

Figure [| shows the typical predictions for hard-scattering cross sections at the Tevatron 
and the LHC, as obtained by using the parton densities of the MRST set. These predic- 
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Gluon Variations That Are Consistent 




Q = 100 GeV 
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Parton x 

Figure 4: A picture of the gluon density uncertainty. The continuous (black) lines refer 
to gluon densities that are constrained only by DIS and DY data. The dashed (coloured) 
lines refer to gluon densities of the MRST group, which uses also prompt-photon data. 

tions have to be supplemented with the corresponding uncertainties |7j] coming from the 
determination of the parton densities and from perturbative corrections beyond the NLO. 
Owing to the increased centre-of-mass energy and to QCD scaling violation (see Fig. |j), 
the kinematic region with small uncertainties is larger at the LHC than at the Tevatron. 

For most of the QCD processes at the LHC, the uncertainty from the parton densities is 
smaller than ±10% and, in particular, it is smaller than the uncertainty from higher-order 
corrections. Some relevant exceptions are the single-jet, W/Z and top quark cross sections. 
In the case of the single-jet inclusive cross section at high {Ex ^ 2 TeV), the uncertainty 
from the poorly known gluon density at high x is larger than that (~ ±10%) from higher- 
order corrections. The W and Z production cross sections are dominated by qq annihilation. 
Since the quark densities are well known, the ensuing uncertainty on the W/Z cross section 
is small (~ ±5%). Nonetheless, in this case the uncertainty from higher-order corrections 
is even smaller, since the partonic cross sections for the DY process are known ]18[ at the 
next-to-next-to-leading order (NNLO) in perturbation theory. In the case of top-quark 
production at the LHC, the gluon channel dominates and leads to an uncertainty of ±10% 
on the total cross section. Also for this process, however, the perturbative component is 
known beyond the NLO. Including all-order resummation of soft-gluon contributions [ |I~9|j , 
the estimated uncertainty from unknown higher-order corrections is approximately ±5% 

00- 
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proton - (anti)proton cross sections 
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Figure 5: QCD predictions for hard-scattering cross sections at the Tevatron and the LHC. 
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3 The gluon density issue 



At present, the processes^ that are, in principle, most sensitive to the gluon density are 
DIS at HERA, 6-quark production at the Tevatron, and prompt-photon production at 
fixed-target experiments. These processes constrain f g for x^lCT 3 , x ~ 10~ 3 -10~ 2 and 
x^lO -1 , respectively. Nonetheless, the gluon density is, in practice, not well determined. 
The issue (or, perhaps, the puzzle) is that from a phenomenological viewpoint the standard 
theory, namely perturbative QCD at NLO, works pretty well for 10~ 3 but not so well 
at larger values of x, while from theoretical arguments we should expect just the opposite 
to happen. This issue is discussed below mainly in its perturbative aspects. We should 
however keep it in mind that all these processes are dominated by hard-scattering scales 
Q of the order of few GeV. Different types of non-perturbative contributions can thus be 
important. 

From the study of DIS at HERA we can extract information on the gluon and sea-quark 
densities of the proton. The main steps in the QCD analysis of the structure functions 
at small values of the Bjorken variable x are the following. The measurement of the 
proton structure function F 2 (x,Q 2 ) ~ Qs{x,Q 2 ) directly determines the sea-quark density 
Is — x (fq + fq)- Then, the DGLAP evolution equation or, more precisely, the following 
equations (the symbol ® denotes the convolution integral with respect to x): 

dF 2 {x,Q 2 )/d In Q 2 ~ P qq ®qs + P qg ®9 , (5) 
dg(x,Q 2 )/d\nQ 2 ~ P gq ®q s + P gg ®9 , (6) 

are used to extract a gluon density g(x,Q 2 ) = xf g (x,Q 2 ) that agrees with the measured 
scaling violation in dF 2 (x, Q 2 ) / d\\iQ 2 (according to Eq. (|5|)) and fulfils the self-consistency 
equation (f|). 

The perturbative-QCD ingredients in this analysis are the AP splitting functions Pabi&s, x). 
Once they are known (and only then), the non-perturbative gluon density can be deter- 
mined. 

The standard perturbative-QCD framework to extract g(x, Q 2 ) consists in using the 
truncation of the AP splitting functions at the NLO. This approach has been extensively 
compared with structure function data over the last few years and it gives a good description 
of the HERA data, down to low values of Q 2 ~ 2 GeV 2 . The NLO QCD fits simply require 
a slightly steep input gluon density at these low momentum scales. Typically IjTTJ , we have 
g(x, Qq) ~ x~ A , with A ~ 0.2 at Qq ~ 2 GeV 2 , and the data constrain g(x, Qq) with an 
uncertainty of approximately ±20%. 

Although it is phenomenologically successful, the NLO approach is not fully satisfactory 
from a theoretical viewpoint. The truncation of the splitting functions at a fixed pertur- 
bative order is equivalent to assuming that the dominant dynamical mechanism leading 
to scaling violations is the evolution of parton cascades with strongly-ordered transverse 
momenta. However, at high energy this evolution takes place over large rapidity intervals 
(Ay ~ lnl/x) and diffusion in transverse momentum becomes relevant. Formally, this 

**The role of jet production at the Tevatron has briefly been recalled in Sect. 0, and it is discussed in 
detail in Ref. (fj. 
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Figure 6: Comparison between the LO (GRV98LO) and NLO (GRV98NLM) GRV 
parametrizations of the gluon and sea-quark densities at Q 2 = 5 GeV 2 . 



implies that higher-order corrections to P a b(as,x) are logarithmically enhanced: 

P ab (a s ,z)~^ + ^ (a s lnx) + ... + — (aslna;)" + ... . (7) 

XX X 

At asymptotically small values of x, resummation of these corrections is mandatory to 
obtain reliable predictions. 



Small- x resummation is, in general, accomplished by the BFKL equation p0|| . In the 
context of structure-function calculations, the BFKL equation provides us with improved 
expressions of the AP splitting functions P a b(as,x), in which the leading logarithmic (LL) 
terms (aslnx)" - , the next-to-leading logarithmic (NLL) terms as{as^x) n , and so forth, 
are systematically summed to all orders n in ag. The present theoretical status of small- x 
resummation is discussed in Ref. |7J. Since in the small- £ region the gluon channel domi- 
nates, only the gluon splitting functions P gg and P gq contain LL contributions. These are 
known p"| pif to be positive but numerically smaller than naively expected (the approach 
to the asymptotic regime is much delayed by cancellations of logarithmic corrections that 
occur at the first perturbative orders in P gg and P gq )- The NLL terms in the quark split- 
ting functions P qg and P qq are known |2^] and turn out to be positive and large. A very 
important progress is the recent calculation |23], Q of the NLL terms in P gg , which are 



found to be negative and large. The complete NLL terms in P gq are still unknown. 

The results of Refs. ||, H, the large size of the NLL terms and the alternating sign 



(from the LL to the NLL order and from the gluon to the quark channel) of the resummed 
small-x contributions have prompted a lot of activity (see the list of references in Ref. 0) on 
the conceptual basis and the phenomenological implications of small- x resummation. This 
activity is still in progress and definite quantitative conclusions on the impact of small- a; 
resummation at HERA cannot be drawn yet. 

At the same time, the capability of the fixed-order approach to produce a good de- 
scription of the proton structure function F 2 (x, Q 2 ) at HERA cannot be used to conclude 
that the small- x behaviour of the gluon density is certainly well determined. In fact, by 
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comparing LO and NLO results, we could argue that the ensuing theoretical uncertainty on 
f g is sizeable [|]]. Going from LO to NLO, we can obtain stable predictions for F 2 , but we 
have to vary the gluon density a lot. As shown in Fig. [|, the NLO gluon density sizeably 
differs from its LO parametrization, not only in absolute normalization but also in x-shape. 
For instance, at x = 10~ 4 and Q 2 = 5 GeV 2 the NLO gluon is a factor of 2 smaller than 
the LO gluon. This can be understood |25[ from the fact that the scaling violation of F 2 
is produced by the convolution P qg ® g (see the right-hand side of Eq. (||)). The quark 
splitting function P qg behaves as 



P qg (a s ,x)~a s P£°)( 



CaOIS 1 
IX X 



(8) 



where the LO term P qg °\x) is flat at small x, whereas the NLO correction is steep. To 
obtain a stable evolution of F 2 , the NLO steepness of P qg has to be compensated by a 
gluon density that is less steep at NLO than at LO. This has to be kept in mind when 
concluding on the importance of small-x resummation because the NLO steepness of P qg 
is the lowest-order manifestation of BFKL dynamics in the quark channel. 

In the large-x region, there is a well-known correlation between as and f g . At small 
x, there is an analogous strong correlation between the x-shapes of P qg and f g . In the 
fixed-order QCD analysis of F 2 , large NLO perturbative corrections at small x can be 
balanced by the extreme flexibility of parton density parametrizations. It is difficult to 
disentangle this correlation between process-dependent perturbative contributions and non- 
perturbative parton densities from the study of a single quantity, as in the case of F 2 at 
HERA. The uncertainty on the gluon density at small x, as estimated from the NLO QCD 
fits of the HERA data, is evidently only a lower limit on the actual uncertainty on f g . 

The production of b quarks at the Tevatron is also sensitive to the gluon density at 
relatively small values of x. The comparison between Tevatron data and perturbative- 
QCD predictions at NLO is shown in Fig. 0. Using standard sets of parton densities, 
the theoretical predictions typically underestimate the measured cross section by a factor of 
2. This certainly is disappointing, although justifiable by the large theoretical uncertainty 
of the perturbative calculation p8 . A lower limit on this uncertainty can be estimated by 



studying the scale dependence and the convergence of the perturbative expansion. Varying 
the factorization and renormalization scales by a factor of four around the 6-quark mass 
m&, the NLO cross section varies by a factor of almost 2 at the Tevatron and by a factor of 
4-5 at the LHC f7j. Similar factors are obtained by considering the ratio of the NLO and 
LO cross sections. 

The present theoretical predictions for 6-quark production at hadron colliders certainly 
need to be improved 0. Since the hard scale Q ~ is not very large, a possible im- 
provement regards estimates of non-perturbative contributions (for instance, effects of the 
fragmentation of the 6-quark and of the intrinsic transverse momentum of the colliding 
partons). As for the evaluation of perturbative contributions at higher orders, the resum- 
mation of logarithmic terms of the type «g \n n (p t /rrib) is important [25] when the transverse 



momentum p t of the b quark is much larger than mj. The resummation of small-x logarith- 
mic contributions «g ln n x can also be relevant, because x ~ 2m^/ y/S is as small as ~ 10~ 3 
at the Tevatron and as ~ 10~ 4 at the LHC. The theoretical tool to perform this resumma- 
tion, namely the fc^-factorization approach p0|, is available. Updated phenomenological 
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Figure 7: Comparison between Tevatron data and NLO QCD for 6-quark production |p7 



The band is obtained by varying factorization and renormalization scales in the NLO 
calculation. 



studies based on this tool and on the information from small-a; DIS at HERA would be 
interesting. 

Prompt-photon production at fixed-target experiments is sensitive to the behaviour 
of the gluon density at large x (x ^0.1). The theoretical predictions for this process, 
however, are not very accurate. Figure ^| shows the factorization- and renormalization- 
scale dependence of the perturbative cross section for the case of the E706 kinematics. If 
the scale is varied by a factor of 4 around the transverse energy of the prompt photon, 
the LO cross section varies by a factor of almost 4. Going to NLO [[Hp the situation 
improves, but not very much, because the NLO cross section still varies by a factor of 
about 2. 



A detailed comparison between NLO QCD calculations and data from the ISR and 
fixed-target experiments has recently been performed in Ref. 



32 1 . As shown in Fig. 



the overall agreement with the theory is not satisfactory, even taking into account the 
uncertainty coming from scale variations in the theoretical predictions. Modifications of 
the gluon density can improve the agreement with some data sets only at the expense of 
having larger disagreement with other data sets. The differences between experiments at 
similar centre-of-mass energies (see, for instance, E706 pBe/530 at v^S* = 31.6 GeV and 
WA70 pp at VS = 23 GeV) are much larger than expected from perturbative scaling 
violations. This can possibly suggest |32| inconsistencies of experimental origin. 



Another (not necessarily alternative) origin of the differences between data and the- 
ory could be the presence of non-perturbative effects that are not included in the NLO 



perturbative calculation. This explanation has been put forward in Refs. |33|, |3J] by in- 
troducing some amount of intrinsic^ transverse momentum (k±) of the colliding partons. 

^To be precise, in Ref. p3 the (kj_) of the colliding partons is not called 'intrinsic', but it is more 
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Figure 8: The dependence on the factorization (/xp) and renormalization (/xr = /x) scale 
of the LO and NLO prompt-photon cross section da/dEx in pN collisions at Et = 5 GeV 
and v^S* = 31.6 GeV. The resummed calculation is discussed in Sect. [|. 



Owing to the steeply falling E T distribution (da/dE T oc l/E T ) of the prompt photon, even 
a small transverse-momentum kick0 can indeed produce a large effect on the cross section, 
in particular, at small values of Et- Phenomenological investigations p3 show that this 



additional (k±) kick can lead to a better agreement between calculations and data. The 
E706 data suggest the value (k±) ~ 1.2 GeV, the WA70 data prefer no (k±), and the UA6 
data in the intermediate range of centre-of-mass energy (VS = 24.3 GeV) may prefer an 
intermediate value of (k±). Similar conclusions are obtained in the analysis by the MRST 
group JIT 



A precise physical understanding of (k±) effects is still missing. On one side, since the 
amount of (k±) suggested by prompt-photon data varies with yS, it is difficult to argue 
that the transverse momentum is really 'intrinsic' and has an entirely non-perturbative 
origin. On the other side, in the case of the inclusive production of a single photon, 
a similar effect cannot be justified by higher-order logarithmic corrections produced by 
perturbative soft-gluon radiation (see Sect. |). A lot of model-dependent assumptions (and 
ensuing uncertainties) certainly enter in the present implementations of the (k±) kick. A 
general framework to consistently include non-perturbative transverse-momentum effects in 
perturbative calculations is not yet available. Recent proposals with this aim are presented 
in Refs. J33| and p6| . 



Further studies on the consistency between different prompt-photon experiments and 



generically called the (fcj_) 'from initial-state soft-gluon radiation'. 

"The Et distribution of the single-photon is not calculable down to Et = or, in other words, dcr/dEr 
is not integrable in the entire kinematic range of Et- Thus, the intrinsic of the incoming partons 
does not simply produce a shift of events from the Iow-Et to the high- Et region. For this reason, the 
terminology l (k±) kick' seems to be more appropriate than l (k±) smearing'. 
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Figure 9: A comparison between NLO QCD calculations and data from the ISR and fixed- 
target experiments for the prompt-photon distribution da/dEx (xt = 2Et/VS). 



on the issue of intrinsic- (k±) effects in hadron-hadron collisions are necessary. Owing to 
the present theoretical (and, possibly, experimental) uncertainties, it is difficult to use 
prompt-photon data to accurately determine the gluon density at large x. Other recent 
theoretical improvements, such as soft-gluon resummation, of the perturbative calculations 
for prompt-photon production at large Xt = are discussed in Sect. [5|. 

Studies of other single-particle inclusive cross sections, such as 7r° cross sections |3 



38], can be valuable to constrain the parton densities and could possibly help to clarify 



some of the experimental and theoretical issues arisen by prompt-photon production. 



4 Partonic cross sections: fixed-order expansions, 
resummed calculations, parton showers 

The calculation of hard-scattering cross sections according to the factorization formula 
(§) requires the knowledge of the partonic cross sections &, besides that of the parton 
densities. The partonic cross sections are usually computed by truncating their perturbative 
expansion at a fixed order in a$: 

*(Pi,P2i Q, {Qi, ■ ■ •}; /4) = «s(/4) {° {LO) (pi, P2; Q, {Qx, ■ ■ •}) (9) 

+ asinl) ( Pl ,p 2 ; Q, {Q u . . .}; /4; /4) 

+ a§(/4) ^ NNLO \p uP2 ; Q, {Q 1 , . . .}; /4 &) + ...}. 

The scale fiR is the arbitrary renormalization scale introduced to define the perturbative 
expansion. Although the 'exact' partonic cross section on the left-hand side of Eq. (|^) 
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does not depend on fi R , each term on the right-hand side (and, hence, any fixed-order 
truncation) separately depends on it. 

The LO (or tree-level) term <r( i0 ) gives only an estimate of the order of magnitude of the 
partonic cross section, because at this order a s is not unambiguously defined. Equivalently, 
we can say that since a^ L °^ does not depend on /i R , the size of its contribution can be varied 
quite arbitrarily by changing fi R in its coefficient a^(fi R ). The strong coupling as can be 
precisely defined only starting from NLO. A 'reliable' estimate of the central value of a thus 
requires the knowledge of (at least) the NLO term a^ NLO \ This term explicitly depends 
on ii R and this dependence begins to compensate that of as([i R ). 

In general, the n-th term in the curly bracket of Eq. (|9]) contains contributions of the 
type (o;s(a*jj) lnQ//iR) n . If fi R is very different from the hard scale Q, these contributions 
become large and spoil the reliability of the truncated expansion (|S|). Thus, in practical 
applications the scale fi R should be set approximately equal to the hard scale Q. As 
mentioned in Sect. |3|, variations of fi R around this central value are typically used to set a 
lower limit on the theoretical uncertainty of the perturbative calculation. 

A better estimate of the accuracy of any perturbative expansion is obtained by con- 
sidering the effect of removing the last perturbative term that has been computed. Since 
«s can be precisely defined only at NLO, this procedure can consistently be applied to 
Eq. (0) only as from its NNLO term. A 'reliable' estimate of the theoretical error on a 
thus requires the knowledge of the NNLO term (ji NNL °) in Eq. (|SJ). 

The LO and NLO approximations of a are used at present in (most of) the fixed-order 
QCD calculations. Prospects towards NNLO calculations of partonic cross sections and 
AP splitting functions are reviewed in Refs. ]7|, |39 . 



The fixed-order expansion (|)|) provides us with a well-defined and systematic frame- 
work to compute the partonic cross section a(pi, p 2 ; Q , {Qi, ■ ■ ■}; fi 2 F ) of any hard-scattering 
process that is sufficiently inclusive or, more precisely, that is defined in an infrared- and 
collinear-safe manner. However, the fixed-order expansion is reliable only when all the kine- 
matical scales Q, {Qi, ■ ■ ■} are of the same order of magnitude. When the hard-scattering 
process involves two (or several) very different scales, say Qi ^> Q, the N n LO term in 
Eq. (|9|) can contain double- and single-logarithmic contributions of the type (asL 2 ) n and 
(asL) n with L = ln(Qi/Q) ^> 1. These terms spoil the reliability of the fixed-order expan- 
sion and have to be summed to all orders by systematically defining logarithmic expansions 
(resummed calculations). 

Typical large logarithms, L = In Q/Qo, are those related to the evolution of the parton 
densities from a low input scale Q to the hard-scattering scale Q. These logarithms 
are produced by collinear radiation from the colliding partons and give single-logarithmic 
contributions. They never explicitly appear in the calculation of the partonic cross section, 
because they are systematically (LO, NLO and so forth) resummed in the evolved parton 
densities f(x,Q 2 ) by using the DGLAP equation (|3]). 

Different large logarithms, L = In Q/\fS ) appear when the centre-of-mass energy V^of 
the collision is much larger than the hard scale Q. These small- x (x = Q/y/S) logarithms 
are produced by multiple radiation over the wide rapidity range that is available at large 
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energy. They usually give single-logarithmic contributions that can be resummed by using 
the BFKL equation. BFKL resummation is relevant to DIS structure functions at small 
values of the Bjorken variable x (see Sect. [^) an d it can also be important at the LHC for 
the production of b quarks and of prompt photons at relatively low Et- 

Another class of large logarithms is associated to the bremsstrahlung spectrum of soft 
gluons. Since soft gluons can be radiated collinearly, they give rise to double-logarithmic 
contributions to the partonic cross section: 

a ~ a**™ jl + £ a n (cW L 2n + C^n-l + C g>_ ^ + ...)} . (10) 

Soft-gluon resummation is discussed in Sect. |^. 

A related approach to evaluate higher-order contributions to the partonic cross sections 
is based on Monte Carlo parton showers (see HOfl and the updated list of references in 



H 0). Rather than computing exactly a^ NLO \ Q-( NNL °) and so forth, the parton shower 
gives an all-order approximation of the partonic cross section in the soft and collinear 
regions. In this respect, the computation of the partonic cross sections performed by 
parton showers is somehow similar to that obtained by soft-gluon resummed calculations. 
There is, however, an important conceptual difference between the two approaches. This 
difference and the limits of applicability of the parton-shower method are briefly recalled 
below. Apart from these limits, parton-shower calculations can give some advantages. 
Multiparton kinematics can be treated exactly. The parton shower can be supplemented 
with models of non-perturbative effects (hadronization, intrinsic k±, soft underlying event) 
to provide a complete description of the hard-scattering process at the hadron level. 

For a given cross section, resummed calculations can in principle be performed to any 
logarithmic accuracy. The logarithmic accuracy achievable by parton showers is instead in- 
trinsically limited by quantum mechanics. The parton-shower algorithms are probabilistic. 
Starting from the LO cross section, the parton shower generates multiparton final states 
according to a probability distribution that approximates the square of the QCD matrix 
elements. The approximation is based on the universal (process-independent) factoriza- 
tion properties of multiparton matrix elements in the soft and collinear limits. Although 
the matrix element does factorize, its square contains quantum interferences, which are 
not positive-definite and, in general, cannot be used to define probability distributions. 
To leading infrared accuracy, this problem is overcome by exploiting QCD coherence (see 
Refs. |ftO| , 4lj and referencees therein): soft gluons radiated at large angle from the 



partons involved in the LO subprocess destructively interfere. This quantum mechanical 
effect can be simply implemented by enforcing an angular-ordering constraint on the phase 
space available for the parton shower evolution. Thus, angular-ordered parton showers can 
consistently compute the first two dominant towers (a§L 2n and oQL 2n ~ l ) of logarithmic 
contributions in Eq. (TO). However, parton showers contain also some subleading logarith- 
mic contributions. For instance, they correctly compute the single-logarithmic terms a^L n 
of purely collinear origin that lead to the LO evolution of the parton densities. Moreover, 
as discussed in Ref. ||3| by a comparison with resummed calculations, in the case of hard- 



scattering processes whose LO subprocess involves two coloured partons (e.g. DIS or DY 
production), angular-ordered parton showers have a higher logarithmic accuracy: they can 
consistently evaluate the LL and NLL terms in Eq. fllBT). The extension of parton-shower 
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algorithms to higher logarithmic accuracy is not necessarily feasible and is, in any case, 
challenging. 



Of course, because of quantum interferences and quantum fluctuations, the probabilistic 
parton-shower approach cannot be used to systematically perform exact calculations at 
NLO, NNLO and so forth. Nonetheless, important progress has been made to include 



matrix element corrections in parton shower algorithms |]44|-|4?|j. The purpose is to consider 
the multiparton configurations generated by parton showering from the LO matrix element 
and to correct them in the hard (non-soft and non-collinear) region by using the exact 
expressions of the higher-order matrix elements. Hard matrix element corrections to parton 
showers have been implemented for top quark decay |38] and for production of W, Z and 



DY lepton pairs [|49], [Kj, [51j| . The same techniques could be applied to other processes, as, 
for instance, production of Higgs boson [^j and vector-boson pairs ]7j . 

Note also that, at present, angular-ordered parton showers cannot be considered as true 
'next-to-leading' tools, even where their logarithmic accuracy is concerned. The consistent 
computation of the first two towers of logarithmic contributions in Eq. (|T0|) is not sufficient 
for this purpose. For instance, to precisely introduce an NLO definition of as, we should 
control all the terms obtained by the replacement as — > as + c a| + 0(ac,). When it is 
introduced in the towers of double-logarithmic terms a^L 2n of Eq. (|TT)P, this replacement 
leads to contributions of the type a^ +1 L 2n ~ a^L 2n ~ 2 . Since these contributions are not 
fully computable at present, the parameter as used in the parton showers corresponds to 
a simple LO parametrization of QCD running coupling. 



5 Soft-gluon resummation 

Double-logarithmic contributions due to soft gluons arise in all the kinematic configurations 
where radiation of real and virtual partons is highly unbalanced (see Ref. |k| and referenes 



therein). For instance, this happens in the case of transverse-momentum distributions at 
low transverse momentum, in the case of hard-scattering production near threshold or when 
the structure of the final state is investigated with high resolution (internal jet structure, 
shape variables). 

Soft-gluon resummation for jet shapes has been extensively studied and applied to 
hadronic final states produced by e + e~ annihilation |l], f|, |54|. Applications to hadron- 
hadron collisions have just begun to appear j55| and have a large, yet uncovered, potential 



(from as determinations to studies of non-perturbative dynamics). 

Transverse-momentum logarithms, L = lnQ 2 /Q 2 ± , occur in the distribution of trans- 
verse momentum Q ± of systems with high mass Q (Q ^> Q±) that are produced with a 
vanishing Q ± in the LO subprocess. Examples of such systems are DY lepton pairs, lepton 
pairs produced by W and Z decay, heavy quark-ant iquark pairs, photon pairs and Higgs 
bosons. In these processes the LO transverse-momentum distribution is sharply peaked 
around Q ± = (da/d 2 Q ± oc ^ 2 >(Q ± )). If the heavy system is produced with Q 2 ± <C Q, 
the emission of real radiation at higher orders is strongly suppressed and cannot balance 
the virtual contributions. The ensuing logarithms, L = lnQ 2 /Q 2 ± , diverge order by order 
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when Q ± — > 0, but after all-order resummation they leads to a finite smearing of the LO 
distribution. 

Threshold logarithms, L = ln(l — x), occur when the tagged final state produced by the 
hard scattering is forced to carry a very large fraction x (x — > 1) of the available centre-of- 
mass energy y/S. Also in this case, the radiative tail of real emission is stronly suppressed 
at higher perturbative orders. Oustanding examples of hard processes near threshold are 
DIS at large x (here x is the Bjorken variable), production of DY lepton pairs with large 
invariant mass Q (x = Q/ \/S), production of heavy quark-antiquark pairs (x = 2rriQ / \/~S) , 
production of single jets and single photons at large transverse energy E T (x = 2Er/VS). 

To emphasize the difference between transverse-momentum logarithms and threshold 
logarithms generated by soft gluons, it can be instructive to consider prompt-photon 
production. In the case of production of a photon paiif\ with invariant mass squared 
Q 2 = {Pi^ +P < 2^) 2 an d total transverse momentum Q ± = p\Q + P2±, transverse-momentum 
logarithms and threshold logarithms appear when Q 2 ± <C Q 2 and Q\ ~ (S/A — Q 2 ), re- 
spectively. However, in the case of production of a single photon with transverse energy 
(or, equivalently, transverse momentum) Et, soft gluons can produce logarithms only in 
the threshold region xt = 2Et/VS — > 1. If the prompt photon has a transverse en- 
ergy that is not closep] to its threshold value, the emission of accompanying radiation is 
not kinematically suppressed and there are no soft logarithms analogous to those in the 
transverse-momentum distribution of a photon pair. In particular, there are no double- 
logarithmic contributions of the type (as In 2 E^/S) n , and perturbative soft gluons are not 
distinguishable from perturbative hard gluons. 

Studies of soft-gluon resummation for transverse-momentum distributions at low trans- 
verse momentum and for hard-scattering production near threshold started two decades 



ago [42, 56]. The physical bases for a systematic all-order summation of the soft-gluon 



contributions are dynamics and kinematics factorizations ||53|| . The first factorization fol- 
lows from gauge invariance and unitarity: in the soft limit, multigluon amplitudes fulfil 
factorization formulae given in terms of universal (process-independent) soft contributions. 
The second factorization regards kinematics and strongly depends on the actual cross sec- 
tion to be evaluated. If, in the appropriate soft limit, the multiparton phase space for this 
cross section can be written in a factorized way, resummation is analytically feasible in 
form of generalized exponentiation of the universal soft contributions that appear in the 
factorization formulae of the QCD amplitudes. 

Note that the phase space depends in a non-trivial way on multigluon configurations 
and, in general, is not factorizable in single-particle contributions^. Moreover, even when 
phase-space factorization is achievable, it does not always occur in the space of the kine- 
matic variables where the cross section is defined. Usually, it is necessary to introduce 
a conjugate space to overcome phase-space constraints. This is the case for transverse- 
momentum distributions and hard-scattering production near threshold. The relevant kine- 



*The same discussion applies to the production of a DY lepton pair. 

^Eventually, when xt <C 1, higher-order corrections are single-logarithmically enhanced. This small- a; 
logarithms, (ashix^)™, have to be taken into account by BFKL resummation. 

■fin the case of jet cross sections, for instance, phase-space factorization depends on the detailed definition 
of jets and it can easily be violated Q. Some jet algorithms, such as the fcj_-algorithm p9| , have better 
factorization properties. 
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matical constraint for (^-distributions is (two-dimensional) transverse-momentum conser- 
vation and it can be factorized by performing a Fourier transformation. Soft-gluon re- 
summation for (^-distributions is thus carried out in 6-space |]60| , |6l| , where the impact 



parameter b is the variable conjugate to Q ± via the Fourier transformation. Analogously, 
the relevant kinematical constraint for hard-scattering production near threshold is (one- 
dimensional) energy conservation and it can be factorized by working in iV-moment space 
|62| , |63| , N being the variable conjugate to the threshold variable x (energy fraction) via a 
Mellin or Laplace transformation. 

Using a short-hand notation, the general structure of the partonic cross section a after 
summation of soft-gluon contributions is 

O = Ores. + <T rem . . (11) 

The term <T res . embodies the all-order resummation, while the remainder <T rem . contains no 
large logarithmic contributions. The latter has the form 

<W = a( f -°-)-[a res .] (f -°- ) , (12) 

and it is obtained from a^'°'\ the truncation of the perturbative expansion for a at a given 
fixed order (LO, NLO, ...), by subtracting the corresponding truncation [<J re s.] ( ' f ' '' > of the 
resummed part. Thus, the expression on the right-hand side of Eq. (|TT|) includes soft-gluon 
logarithms to all orders and it is matched to the exact (with no logarithmic approximation) 
fixed-order calculation. It represents an improved perturbative calculation that is every- 
where as good as the fixed-order result, and much better in the kinematics regions where the 
soft-gluon logarithms become large (ctsL ~ 1). Eventually, when agL^l, the resummed 
perturbative contributions are of the same size as the non-perturbative contributions and 
the effect of the latter has to be implemented in the resummed calculation. 

The resummed cross section has the following typical form: 

a res .=a k s [ &M.C-S, (13) 

J inv. 

where the integral J inv denotes the inverse tranformation from the conjugate space where 
resummation is actually carried out. Methods to perform the inverse transformation are 
discussed in Refs. and for Q±-resummation and threshold resummation, respec- 



tively. The C term has the perturbative expansion 

C = 1 + d a s + C 2 a 2 s + . . . (14) 

and contains all the constant contributions in the limit L — > oo (the coefficients Ci, C%, ■ ■ ■ 
do not depend on the conjugate variable). The singular dependence on L (more precisely, 
on the logarithm L of the conjugate variable) is entirely exponentiated in the factor S: 

S = exp {L gi(a$L) + g 2 (a s L) + a s g 3 (a s L) + . . .} . (15) 

In the exponent, the function Lgi resums all the leading logarithmic (LL) contributions 
agL n+1 , while g 2 contains the next-to-leading logarithmic (NLL) terms a§L n and so forth| 

§To compare this notation with that of Ref. fl68fl , we can notice that our functions en a re obtained by 
the straightforward integration over ~p, of the functions A(as(jT)) and B(as('p)) of Ref. |p8| . In particular, 
our terms 31, (72, ff3 are not to be confused with the non-perturbative parameters of the same name used in 
Ref. M. 
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(all the functions are normalized as g>j(A = 0) = 0). Note that the LL terms are formally 
suppressed by a power of as with respect to the NLL terms, and so forth for the successive 
classes of logarithmic terms. Thus, this logarithmic expansion is as systematic as the fixed- 
order expansion in Eq. (|9]). In particular, using a matched NLL+NLO calculation, we 
can consistently i) introduce a precise definition (say MS) of «s(a*) and ii) investigate the 
theoretical accuracy of the calculation by studying its renormalization-scale dependence. 

The structure of the exponentiated resummed calculations discussed so far has to be 
contrasted with that obtained by organizing the logarithmic expansion on the right-hand 
side of Eq. (|l^) in terms of towers as 

a ~ a k s a^ {h^L 2 ) + a s L t 2 (a s L 2 ) + a 2 s L 2 t 3 (a s L 2 ) + . . .} , (16) 

where the double-logarithmic function t^a^L 2 ) and the successive functions are normalized 
as U(0) = const. While the ratio of two successive terms in the exponent of Eq. ( |15D is 
formally of the order of as, the ratio of two successive towers in Eq. ( fl6|) is formally of 
the order of a$L. In other words, the tower expansion sums the double-logarithmic terms 
(a>sL 2 ) n , then the terms «gL 2n_1 ~ asL(asL 2 ) n ~ l , and so forth; it thus assumes that the 
resummation procedure is carried out with respect to the large parameter a$L 2 (a^L 2 ^ 1). 
On the contrary, in Eq. (|15|) the large parameter is a$L&l. The tower expansion allows 
us to formally extend the applicability of perturbative QCD to the region L ^ l/y/a s , and 
exponentiation extends it to the wider region L^l/ag. This fact can also be argued by 
comparing the amount of information on the logarithmic terms that is included in the 
truncation of Eqs. (|Uj) and ([16]) at some logarithmic accuracy. The reader can easily check 
that, after matching to the NLO (LO) calculation as in Eq. (|TTD , the NLL (LL) result of 
Eq. (fT5|) contains all the logarithms of the first four (two) towers in Eq. flT6|) (and many 
more logarithmic terms). 

In the case of Q_|_-distributions, full NLL resummation has been performed for lepton 
pairs, W and Z bosons produced by the DY mechanism |66| and for Higgs bosons 



produced by gluon fusion JB7J. Corresponding resummed calculations are discussed in 
Refs. |52| , IBlj and references therein. 

Threshold logarithms in hadron collisions have been resummed to NLL accuracy for DIS 
and DY production |62], ^H, [H, |69j and for Higgs boson production |7(| . Recent theoretical 



progress [[TT], [F2|, [19| regards the extension of NLL resummation to processes produced by LO 
hard-scattering of more than two coloured partons, such as heavy-quark hadroproduction 
7T| , [0| and leptoproduction |75| , as well as prompt-photon |7^-|77|, quarkonium JTB] and 



vector-boson [T9j production. 

An important feature of threshold resummation is that the resummed soft-gluon con- 
tributions regard the partonic cross section rather than the hadronic cross section. This 
fact has two main consequences: i) soft-gluon contributions can be sizeable long before the 
threshold region in the hadronic cross section is actually approached, and ii) the resum- 
mation effects typically enhance the fixed-order perturbative calculations. 

The first consequence follows from the fact that the evolution of the parton densities 
sizeably reduces the energy that is available in the partonic hard-scattering subprocess. 
Thus, the partonic cross section a in the factorization formula (§) is typically evaluated 
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much closer to threshold than the hadronic cross section. In other words, the parton 
densities are strongly suppressed at large x (typically, when x — > 1, f(x, /i 2 ) ~ (1 — x) v 
with 1] ~ 3 and rj ~ 6 for valence quarks and sea-quarks or gluons, respectively); after 
integration over them, the dominant value of the square of the partonic centre-of-mass 
energy (s) = (x\X2)S is therefore substantially smaller than the corresponding hadronic 
value S. 

The second consequence, which depends on the actual definition of the parton densities, 
follows from the fact that the resummed contributions are those soft-gluon effects that are 
left at the partonic level after factorization of the parton densities. After having absorbed 
part of the full soft-gluon contributions in the customary definitions (for instance, those 
in the MS or DIS factorization schemes) of the parton densities, it turns out that the 
residual effect in the partonic cross section is positive and tends to enhance the perturbative 
predictions. 

A quantitative illustration of these consequences is given below by discussing top-quark 
and prompt-photon production. The discussion also shows another relevant feature of 
NLO+NLL calculations, namely, their increased stability with respect to scale variations. 

The effects of soft-gluon resummation on the top-quark production cross sections at 



hadron colliders have been studied in Refs. [19, pQ-p3fl. In the case of pp collisions, the 



comparison between QCD predictions at NLO and those after NLL resummation is shown 
in Fig. W |BJ. At the Tevatron the resummation effects are not very large and the NLO 



cross section is increased by about 4%. This had to be expected because the top quark is not 
produced very close to threshold (x = 2m t /VS ~ 0.2, at the Tevatron). Note, however, that 
the dependence on the factorization/renormalization scale of the theoretical cross section 
is reduced by a factor of almost 2 by including NLL resummation. More precisely, the scale 
dependence (~ ±5%) of the NLO+NLL calculation becomes comparable to that obtained 



by using different sets of parton densities |TTJ . Combining linearly scale and parton density 



uncertainties, the NLO+NLL cross section is o t i = 5.0 ± 0.6, with m t = 175 GeV and 
VS = 1.8 TeV ||. 

At the LHC (x = 2m t / \/S ~ 0.03) the top quark is produced less close to the hadronic 
threshold than at the Tevatron. However this is compensated by the fact that the gluon 
channel^ is more important at the LHC. As a result, the effect of including soft-gluon 
resummation to NLL accuracy is very similar: the NLO cross section is enhanced by ~ 5% 
and its scale dependence is reduced from ~ ±10% to ~ ±5%. Note, however, that the 
uncertainty (~ ±10%) coming from the parton (gluon) densities is larger than at the 
Tevatron [[7]. 



Similar qualitative results are obtained |76j when NLL resummation is applied to 



prompt-photon production at fixed-target experiments. The scale dependence of the the- 
oretical calculation is highly reduced and the resummed NLL contributions lead to large 
corrections at high xt = 2Et/vS (and smaller corrections at lower xt)- Of course, the 
impact of soft-gluon resummation is quantitatively more sizeable in prompt-photon pro- 
vince f g is steeper than f q at large x, partonic cross sections in gluon subprocesses are typically closer 
to threshold than in quark subprocesses. Moreover, the intensity of soft-gluon radiation from gluons is 
larger than that from quarks by a factor of ~ Ca/Cf ~ 2. 
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Figure 10: The ti production cross section in pp collisions as a function of v^S*- The 
solid lines represent the NLO results for different choices (fi = m t /2 and /x = 2m t ) of the 
renormalization/factorization scale /i — /ir — (j,f, normalized to the result with \i = m t . 
The dashed lines represent the NLO+NLL results for different choices of /i (/i = m t /2,m t 
and 2m t ), normalized to the NLO result with \i = m t . 
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Figure 11: The dependence on the factorization {hf) and renormalization = /i) scale 
of the prompt-photon cross section da/dE? in pN collisions at Ex = 10 GeV and = 
31.6 GeV. The short-dashed, long-dashed and solid lines are respectively the results at LO, 
NLO and after NLO+NLL resummation. 
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Figure 12: The prompt-photon cross section da/dEx in pN collisions at Ex = 10 GeV 
as a function of the energy -Ebeam of the proton beam. The solid lines represent the NLO 
results for different choices (/x = Ex/2 and /x = 2Ex) of the renormalization/factorization 
scale /x = fiR = up, normalized to the result with /x = Ex- The dashed lines represent the 
NLO+NLL results for different choices of /x (xx — Ex/2, Ex and 2Ex), normalized to the 
NLO result with /x = Ex- 



duction than in top-quark production, because xx can be as large as 0.6, the hard scale 
Ex is much smaller than m t (thus, as(Ex) > as{m t )) and the gluon channel is always 
important. The scale dependence of the theoretical cross section for the E706 kinematics 



is shown in Fig. |TT]. Fixing = /j,p = /j and varying fi in the range Ex/2 < /x < 2 Ex with 
Ex = 10 GeV, the cross section varies by a factor of ~ 6 at LO, by a factor of ~ 4 at NLO 
and by a factor of ~ 1.3 after NLL resummation. The highly reduced scale dependence 



of the NLO+NLL cross section is also visible in Fig. |H| which, in particular, shows that 
when Ex = 10 GeV and E^ ea _ m = 530 GeV the central value (i.e. with /x = Ex) of the NLO 
cross section increases by a factor of ~ 2.5 after NLL resummation. As expected, the size 
of these effects is reduced by increasing y/S at fixed E T (see Fig. |l^) or by decreasing E T 
at fixed \fS (see Fig. |]). 

The comparison with the E706 data shown in Fig. [13] suggests that the NLO+NLL 
calculation can help to better understand prompt-photon production at large xx- Note, 



however, that this comparison has to be regarded as preliminary in several respects [76|. In 



particular, the parton densities used in Fig. O are those extracted from NLO fits. Owing 



to the soft-gluon enhancement at large Xx, refitting the parton densities may lead to a 
smaller f g at large x and, consequently (because of the momentum sum rule), a larger f g 
at intermediate x. As a result, this procedure could somehow increase the theoretical cross 
section also at smaller values of xx- 

Soft-gluon resummation at NLL accuracy is now available for all the processes (namely, 
DIS, DY and prompt-photon production) that are typically used to perform global fits to 
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Figure 13: E706 prompt-photon data compared with theoretical calculations, which use the 
parton densities of the set CTEQ4 and GRV photon fragmentation functions. The solid 
and dashed lines correspond to the NLO+NLL and pure NLO calculations, respectively. 

parton densities. A detailed extraction/evolution of parton densities by consistently using 
NLL resummed calculations is thus nowadays feasible. 



6 Other topics 

The activity of the QCD Working Group at this Workshop has also been devoted to other 
topics, such as automatic computation of matrix elements and LO cross sections for mul- 
tiparticle processes at high-energy colliders, definition and properties of jets algorithms, 
definition of isolated photons and related NLO calculations. Corresponding contributions 
are included in these Proceedings. 

Other studies performed during this Workshop have a large overlap with the activity of 
the related Workshops at FERMILAB and CERN and can be found in those Proceedings 

Acknowledgements. I would like to thank the members of the Local Organizing Com- 
mittee for the excellent Workshop. 
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